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\Q . In terms of the heavy chiral Lagrangian and the unitarized coupled-channel scat- 

o ■ 

I tering amplitude, interaction between the heavy meson and the light pseudoscalar 

meson is studied. By looking for the pole of scattering matrix on an appropriate 
1) . Riemann sheet, a DK bound state D*q with the mass of 2.312 it 0.041 GeV is found. 

^ I This state can be associated as the narrow L'*j(2317) state found recently. In the 

^ ' same way, a BK bound state B*q is found, and its mass of 5.725 it 0.039 GeV is 

H : 

- - ■ predicted. The spectra of Dq and Bq with 7=1/2 are further investigated. One 

broad and one narrow states are predicted in both charm and bottom sectors. The 
coupling constants and decay widths of the predicted states are also calculated. 
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I. INTRODUCTION 



The recently discovered narrow- width state P*t(2317) 1| stimulates both experimental 



and theoretical 




ill 



that this new state is a conventional cs state 





interest. Many physicists surmised 

y, H H Q Q Q, Q, Q 



and the others believed that it can be an exotic meson state, such as a four-quar 



19,0, 21, 22,0 



c state 



a Df,7T quasi-bound state f24|, a DK bound state 
mixed state of cs with DK |29l. l3n| or with four-quark state 0,0], and etc. On the other 
hand, one proposed that 7(2317) with = 0+ could be the chiral partner of the ground 

nn n 

state of Ds P, 11^- However, the author in Ref. mentioned that the chiral doubler 
produced by using Random Phase Approximation equations should be (Ds(1968),Ds(2392)) 
rather than {Ds{1968),Ds{23n)), although the scalar state Z)*(2392), as the scalar chiral 
partner of Ds(1968) state, has not been found yet I33I. Up to now, the structure of D*j(2317) 
is still indistinct and should carefully be studied. Moreover, the Belle collaboration recently 
reported a broad 0^ charmed meson with mass and width being m^*o = 2308 ± 60 MeV and 
r^.o = 276 ± 99 MeV, respectively 0|, and the FOCUS collaboration reported a broad 0"*" 
charmed meson with mass and width being mn.o = 2407 ± 56 MeV and Fnto = 240 ±114 

. Though they are consistent with each other within experimental 



errors, whether they are the same particle is still in dispute 



On the other hand, it has been shown that the light scalar mesons a, /o(980), ao(980) and 
K can dynamically be generated through the .S* wave interaction between Goldstone bosons 

45! |4^. In such an approach. 



in the chiral unitary approach (ChUA) 0,0,00 

the amplitudes from the chiral perturbation theory (ChPT) are usually adopted as the 
kernels of the factorized coupled- channel Bethe-Salpeter (BS) equations. In this procedure, a 
Lagrangian in a specific expanded order, where the symmetries of ChPT should be preserved, 
is chosen at the beginning, and then the higher order corrections to the amplitudes are re- 
summed with the symmetries kept up to the order of the expansion considered. Namely, 
what the unitary CHPT does in the successive step is re-summing a string of infinite loop 



47 



m 



49|. Moreover, ChUA has 



diagrams while the the symmetries of ChPT are held 
been applied to study the S wave interaction between the lower lying vector meson and 
the Goldstone boson, and most of the known axial-vector mesons can also be generated 



dynamically js^. Based on the valuable achievements mentioned above, extending ChUA to 
the heavy-light meson sector to study the 5* wave interaction between the heavy pseudoscalar 
meson and the Goldstone boson, and consequently the structures of possible heavy scalar 
mesons, would be extremely meaningful. In fact, similar work, called x-BS(3) approach, 
has been done [2^, 22|. In such an approach, heavy-light meson resonances and open-charm 
meson resonances were predicted through checking speed plots together with the real and 
imaginary parts of the reduced scattering amplitudes. In our opinion, studying the poles on 
the appropriate Riemann sheet of the scattering amplitude would be a powerful procedure 
to reveal the properties of the generated states in a more accurate way. In this paper, the 
S wave interaction between the heavy meson and the light pseudoscalar meson is studied 
by using the extended chiral unitary approach, called heavy chiral unitary approach. The 
poles that associate with the experimentally observed narrow D*j{2317) and broad Dq in 
the / = 0, S = 1 and I = S = channels, where / and S denote the isospin and the 
strangeness, respectively, are searched. The corresponding coupling constants and decay 
widths are also discussed. 



II. COUPLED-CHANNEL HEAVY CHIRAL UNITARY APPROACH 



In order to describe the interaction between the Goldstone boson and the heavy pseu- 
doscalar boson, we employ a leading order heavy chiral Lagrangian 0, |^ 

1 



J TT 



(1) 



where /^r = 92.4 MeV is the pion decay constant, P represents the charmed mesons 
(cm, cd, cs), namely (Z)°, D^, Df), and $ denotes the octet Goldstone bosons and can be 
written in the form of 3 x 3 matrix 



V 



TT 

K- 



V2" ^ V6' 



(2) 



This Lagrangian is equivalent to the SU (4) extrapolation of the ordinary meson meson chiral 
Lagrangian, eliminating the exchanges of heavy vector mesons in the equivalent picture of 
vector meson exchange [54]. Obviously, the similar investigation in the bottom sector can 
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be carried out by replacing P in Eq. (^J with the anti-bottom mesons {bu, bd, bs), namely 

We are interested in the heavy mesons in the / = 0, S = 1 and I = ^, S = channels 
that can be specified by their own isospins, respectively. In terms of Eq. (Q), the amplitudes 
can easily be obtained by 

Vi{s,t,u) = ^^{s-u), (3) 

where i and j represent the initial state and the final state, respectively. In the / = case, 
i (j) can be 1 and 2 which represent the coupled DK and Dgf] channels in the charmed 
sector, respectively, and BK and Bgrj channels in the bottom sector, respectively. In the 
J = I case, i (j) can take 1, 2 and 3 which denote the coupled Dtt, Drj and DgK channels 
in the charmed sector, respectively, and Bn, Brj and BgK channels in the bottom sector, 
respectively. The coefficients C/ are listed in Table HI 



TABLE I: Coefficients C/- in Eq. ©. 



r^O /^O r^O 
^11 "-"12 "-"22 


^1/2 ^1/2 ^1/2 ^1/2 
"-'11 "-"12 "-"22 "-'13 


"-'23 


^1/2 
"-"33 


-2^/3 


-2 


2 


-1 



The tree level amplitudes can be projected to the S wave by using 

1 



V/j' = - J ^dcosev/j{s,t{s,cose),u{s,cose)), (4) 



and 



2 2 o /r 2 , Ks,rnl,ml) \{s,mlml 
s — 7712 — mi — 2\/ [mi H J [m^ + 



As 4s 
+ — \/A(s,mf,mi)A(s,m§,m|) cos 6*, (5) 

where X{s,m'^,m^) = [s — {rrii + mj)'^][s — {rrii — rrij)'^] and the on-shell condition for the 
Mandelstam variables, s + t + u = Yli=i f^fi is applied. 

In ChUA, under the on-shell approximation, the full scattering amplitude can be con- 
verted into an algebraic BS equation 

T = (1 - VG)- V, (6) 
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where V is a matrix whose elements are the 5* wave projections of the tree diagram amph- 
tudes and G is a diagonal matrix with the element being a two-meson loop integral 

where pi and p2 are the four-momenta of the two initial particles, respectively, and mi and 



m2 are the masses of the particles appearing in the loop. It was s^ 



lown that the scattering 



matrix derived in such a way satisfies the unitary relation ^,42, 4fi| . 

The loop integral can usually be calculated in the center-of-mass frame by using a three- 
momentum cut-off parameter qmax 13 • However, in this method, an artificial singularity 
of the loop function might be produced and the applicability of the method is limited. 
The better way to remove the singularity of the loop integral is using the dispersion relation 
where a subtraction constant is employed. Then, the analytic expression of Gu^s) can be 
expressed by j42| 

+ ^[log (s - A + (t) + log (s + A + a) 

-log(-s + A + (t) -log(-s- A + o-)]}, (8) 

where a(/i) is the subtraction constant, fi denotes the regularization scale, a = [— (s — (mi + 
^2Y){s — (mi — m2Y)Y^'^ and A = m\ — m^. This result is independent of /i, because the 
change in Ga, caused by a variation of /i, is cancelled by the corresponding change of the 
subtraction constant a(/i). 

III. POLES ON APPROPRIATE RIEMANN SHEETS 

The physical states are closely associated with the poles of the scattering amplitude on the 
appropriate Riemann sheet of the energy plane. For instance, considering only one channel, 
a bound state is associated with a pole below the threshold value in the real axis of the 
energy plane, and the three-momentum of the scattered meson in the center of mass frame 
of the two mesons system can be written as Pcm = AVcml- A resonance should be related 
with a pole on the second Riemann sheet, namely, Impcm < 0. In the coupled channel case, 
the situation is somewhat complicated. Detailed relation can be found in Ref. 



6 



Before searching for poles of the scattering amphtude, the range of subtraction constant 
values in the dispersion relation method should firstly be estimated. It can be done by 
comparing the calculated value of loop integration in the dispersion relation method with 
the one obtained in the cut-off method, although there might be an artificial singularity 
problem in the cut-off method 46]. The cut-off momentum can approximately be chosen as 

qmax ~ -m^, (9) 

where is the mass of the Goldstone boson and denotes the chiral symmetry breaking 
scale which is about 1 GeV. The resultant Qmax for (p = it, K and t] are all in the region of 
0.8-0.9 GeV. Thus, it is reasonable to pick up a value of qmax in the region of 0.8 ± 0.2 GeV. 
Then, we adjust the renormalization scale /i or the subtraction constant a(yu) to match 
the calculated value of the loop integral in the dispersion relation method with the one 
obtained in the cut-off method at y/s = m£,{TnB) + ttik in a specific qmax value case, say 
qmax = 0.6, 0.8 and 1.0 GeV, respectively. The resultant loop integration curves versus s in 
two different methods are very close in the region around and below the matching point ^/s. 
The corresponding values of a(/i) and qmax are tabulated in Table |Tl| With the estimated 
a(/i) value, the full scattering amplitude can be calculated. 

TABLE II: The values of a(^) from matching. We use n = mo for the charm sector, and ^ = ms 
for the bottom sector, respectively. 



Qmax (GeV) 


0.6 


0.8 


1.0 


o("t-d) 


-0.373 


-0.630 


-0.864 


a{mB) 


0.0232 


-0.0856 


-0.187 



The poles of the scattering matrix in the J = 0, S = 1 channel in both the charmed 
sector and bottom sector are searched for first. It is shown that on the first Riemann sheet 
of the energy plane, there is only one pole located on the real axis below the lowest strong 
decay threshold, m£, + itik = 2.367 GeV, in the charmed sector and only one pole on the 
real axis below the lowest strong decay threshold, uib + ttlk = 5.773 GeV, in the bottom 
sector as well. The resultant pole positions with different which correspond to the 

qmax = 0.6, 0.8 and 1.0 GeV cases, are tabulated in Table IIIIl respectively. These poles 
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TABLE III: Poles in the (/, S) = (0, 1) channel. 



Qmax 


(GeV) 


0.6 


0.8 


1.0 




(GeV) 


2.353 


2.317 


2.270 


BsO 


(GeV) 


5.764 


5.729 


5.661 



are apparently associated with the DK bound state and the BK bound state, respectively. 
Due to the existence of the s quark, these bound states should be scalar heavy mesons, 
namely D*q and B*q, respectively. More specifically, when a{m£,) = —0.630, corresponding 
to qmax = 0.8 GeV, the mass of the DK state, namely D*^ , is about 2317 MeV, which is 
almost the same as the measured value of D*j(2317). Taking into account the uncertainty 
of subtraction constant, the mass of the D*q (0, 1) state in our model is 2.312 ± 0.041 GeV. 
Also due to the uncertainty of 0(771^), the predicted mass of the BK bound state, namely 
Btn (0, 1) state, is 5.725 ± 0.039 GeV. This mass is consistent with the mass predicted in 

Ti n n n 

Refs. 121123], but larger than that in Refs. |la,|2a]. For comparison, we list the mass of B*^ 
predicted in different models in Table HVl 



TABLE IV: Mass of B*q predicted in different models. 





Our result 






m 


m 


ms:, (GeV) 


5.725 ± 0.039 


5.728 ±0.035 


5.71 ±0.03 


5.627 


5.643 



In the / = |, S = case, the poles are located on nonphysical Riemann sheets. Usually, 
if Impcm is negative for all the channels open for a certain energy, the width obtained would 
correspond more closely with the physical one. We search for poles in this particular sheet. 

There are two poles in either charmed sector or bottom sector. The width of the lower pole 
is broad and the width of the higher one is narrow. The obtained poles are listed in Table 
IVI In either the charmed or bottom sector, the lower pole is located on the second Riemann 
sheet {Impcmi < 0, lmpcm2 > 0, Impcms > 0, where Pcmi denotes the momentum of one of 
the interacting mesons in the i-th channel in the center of mass system). This pole should 
be associated with a Dtt {Bit) resonance in the charmed (bottom) sector. Consequently, 
this state should easily decay into Dti {Bn) in the charmed (bottom) sector. 
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TABLE V: Poles in (/, 5) = (i, 0) channel. 



qmax (GeV) 


0.6 


0.^ 




1.0 






2.115 


- f0.147 


2.099 - 


fO.lOO 


2.079 - 


f0.067 


Dl (GeV) 














2.488 


- iiO.039 


2.445 - 


i0.049 


2.429 - 


i0.002 


So* (GeV) 


5.564 
5.864 


-«0.160 
- i0.027 


5.534 - 
5.827 - 


iO.llO 
i0.026 


5.507 - 
5.821 - 


i0.074 
i0.019 



The higher pole in either charmed or bottom sector is found on the third Riemann sheet 
(Impcmi < 0, Impcm2 < 0, Impcms > 0) when a(yu) corresponds to q^ax = 0.6 GeV or 0.8 
GeV, or on the second Riemann sheet when a(/i) corresponds to qmax = 1.0 GeV. The pole 
should be associated with an unstable DsK {BgK) bound state in the charmed (bottom) 
sector due to its narrow width. It should be mentioned that the situation for the higher pole 
in the later case, namely a(/i) corresponding to qmax = 1.0 GeV, is somewhat complicated. 
Besides a pole on the second Riemann sheet, poleji = 2.429 — z0.002 GeV shown in Table 
El there is a shadow pole, polem = 2.397 — ^0.043 GeV, on the third Riemann sheet. 
Note that Re(po/e//) > rriD + rrir, and Re(po/e///) < mo + rrirj. A sketch plot for the 
paths of these two poles to the physical region in the energy plane is shown in Fig. 
From this cartoon, one sees that poleu corresponds more closely with the physical one. 
Therefore, we choose polejj = 2.429 — z0.002 GeV as the result. Similar complexity appears 
at polem = 2.488 - i0.039 GeV in Table El due to the existence of poley = 2.048 - i0.020 
GeV. With the same reason, we disregard poley- 

Physical region 

FIG. 1: Paths from polejj on Riemann sheet II and polem on Riemann sheet III to the physical 
region in the energy plane, where Ei = mo + and E2 = mo + mrj. 



Considering the deviations of the data caused by the uncertainty of a(/x) Table El we 
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predict the mass and the width of the broad (i, 0) state as 2.097 ± 0.018 GeV and 
0.213 ± 0.080 GeV, respectively, and the mass and the width of the narrow D* (|, 0) state 
as 2.448 ± 0.030 GeV and 0.051 ± 0.047 GeV, respectively. In the same way, we forecast the 
mass and the width of the broad B* (i, 0) state as 5.536 ± 0.029 GeV and 0.234 ± 0.086 
GeV, respectively, and the mass and the width of the narrow B* (i, 0) state as 5.842 ±0.022 
GeV and 0.035 ± 0.019 GeV, respective 
Recalling the predictions in Refs. 



m |27|, we noticed that by checking the reduced 



scattering amplitude curves in the speed plot, the authors in Ref. j26| found a broad 
state with mass of 2138 MeV and a narrow states with mass of 2413 MeV in the charmed 
sector, and by further adjusting free parameters in the next-to-leading order to reproduce 
the i:'*o(2317) state with mass of 2317±3 MeV and the D* state with mass of 2308±60 MeV 



and width of 276 ±99 MeV given in Ref. [37[, the authors in Ref. j27l] obtained a broad state 
with mass of 2255 MeV and width of about 360 MeV and predicted a very narrow state with 
mass of 2389 MeV. In the same way, the authors in Ref. j2a| further predicted a broad state 
with mass of 5526 MeV and a narrow states with mass of 5760 MeV and width of about 30 
MeV in the bottom sector. It seems that our predicted Dq (|, 0) states are consistent with 



37 



It should 



those in Ref. j27|], although they still deviate from the experimental data 
be mentioned that because of the large uncertainty in the data analysis and existence of the 
predicted higher narrow state just around the D2(2460) region, the present model could not 
be disregarded rudely. 



IV. COUPLING CONSTANTS AND DECAY WIDTHS 



The decay properties of predicted states are studied by making the Laurent expansion of 
the amplitude around the pole 5^ 

rr _ 9i9j 



± 70 ± ll{S - Spole) + 



(10) 



where Qi and gj are coupling constants of the generated state to the i-th and j-th channels. 
giQj can be obtained by calculating the residue of the pole 



gigj = lim (s - Spoie)Tij . 

S ^ Spole 



(11) 



In the case where a(/i) corresponds to qmax = 0.8 GeV, we calculate the residues of the 
poles, and consequently the coupling constants. The resultant coupling constants for the 
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D*o and B*Q {D^ and B^) states are tabulated in Table EH ()VII|) . From these tables, one 
sees that the coupling constants again are consistent with the results in the pole analysis. 
In the (0, 1) channel, the coupling of D*q {B*^) to the DsT] {Bs1]) channel is weaker than 
that to the DK (BK) channel. This is because the D*q {B*q) state is the DK (BK) bound 
state. In the (|, 0) channel, the coupling of the lower broad (Sq) state to the Dn (Bn) 
channel is stronger than the coupling of the higher narrow one to the Dtt (Btt) channel; the 
coupling of the lower state to the Dn {Bn) channel is stronger than that to the DgK {BgK) 
channel and the Dt] {Brf) channel, and the coupling of the higher state to the D^K [BgK) 
channel is stronger than that to the Dr] [Br]) channel and the Dn {Bn) channel. These 
are consistent with the pole analysis for the lower pole being a Dn (Bit) resonance and the 
higher pole being the unstable bound state of DgK {BgK). 

TABLE VI: Coupling constants of the generated D*q and B*q states to relevant coupled channels. 
In this case, gi and g2 are real. All units are in GeV. 





Masses 


bil 


152 1 




2.317 


10.203 


5.876 




5.729 


23.442 


13.308 



TABLE VII: Coupling constants of the generated and states to relevant coupled channels. 
All units are in GeV. 





Poles 


91 


\9i\ 


92 


I52I 


53 I53I 




2.099 


- iO.lOO 


7.750 + ^5.191 


9.328 


-0.184 + i0.096 


0.208 


4.648 + i3.083 5.578 


Dl 2.445 


- i0.049 


0.030 + i3.636 


3.636 


-6.845 - i2.248 


7.205 


-10.815 + il. 543 10.924 


Bl 


5.534 


- iO.llO 21.443 + Z12.060 24.602 


-2.239 - f0.730 


2.355 


13.503 + «7.016 15.217 




5.827 


- i0.026 


0.256 + ^6.958 


6.963 


-14.697 -i4.880 15.486 


-25.000 - i0.602 25.003 



The decay widths of generated states are further evaluated. We first study the states in 
the (0, 1) channel. The D*q state cannot decay into either DK or D^t], because the mass of 
the state is lower than the threshold of the DK channel. Moreover, the D*q'(2317) —>■ D^-k^ 
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decay violates the isospin symmetry. Thus, the decay width of D*o'(2317) should be very 
small. This decay can only occur through n^-r] mixing. According to Dashen's theorem jSSf , 
the TT^-r] transition matrix should be 

t^^ = {7r^\n\r]) = -0.003 GeV, (12) 

and the decay width reads 



Pcm I 5'2^7rr) |2 



(13) 



V 



where M is the mass of the initial state, g2 represents the coupling of D*q(2317) to Dt], and 
Pcm denotes the three-momentum in the center of mass frame and can be written as 

Pcm = ^ ^{M^ - (mz)+ + m^o)2)(M2 - {mo+ - m^o)2). (14) 

Then, the partial decay width of the D*(J"(2317) D^t^^ process can be obtained as 

r(D*+(2317) ^ D+7r°) = 8.69 keV. (15) 

This value is compatible with that in Ref. Similarly, the partial decay width of the 

isospin violated decay i?*Q(5729) — > -B°7r° can be evaluated as 

r(5;°(5729) ^ S°7r°) = 7.92 keV. (16) 

We then study the states in the (|, 0) channel. For the higher state, two strong decay 
channels are opened. The fraction ratio of the decay widths for these two decay channels 
can be calculated by utilizing the coupling constants given in Table I VI II Let Fi and F2 
denote the partial decay widths with the final states being D{B)7i and D{B)r], respectively. 
The ratio Fi/(Fi + F2) can be written by 



For higher Dq and Bq states, we have 

R{D*q) = 0.446, R{B*) = 0.829. (18) 

It is shown that in the bottom sector, the higher narrow state is easier to decay into Bn 
than into Brj, but in the charmed sector, the higher narrow state can decay into Dir and 
Dt] in almost the same weight. 
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V. CONCLUSION 



Based on the heavy chiral unitary approach, the S wave interaction between the pseu- 
doscalar heavy meson and the Goldstone boson is studied. By calculating full scattering 
amplitudes via an algebraic BS equation, the poles on some appropriate Riemann sheets are 
found. These poles can be associated with bound states or resonances. With a reasonably 
estimated single parameter a(/i) in the loop integration, a pole on the real axis on the first 
Riemann sheet, which is associated with the bound state, in the two-coupled-channel calcu- 
lation in the (0, 1) channel is found. Because the mass of the pole in the charmed sector is 
about 2.312 ± 0.041 GeV, this state should be a 0+ DK bound state and can be regarded 
as the recently observed D*j(2317). Meanwhile, a 0+ state B*q, which should be a BK 
bound state, is predicted. Its mass is about 5.725 ± 0.039 GeV. In the I = \-, S = case, 
three-coupled-channel calculations are performed in both charmed and bottom sectors. In 
the charm sector, a broad pole structure, which is associated with a resonance, is found at 
about (2.097 ± 0.018 - i0.107± 0.040) GeV. Besides, a narrow pole structure, which can be 
interpreted as a quasi-bound state of D^K, at about (2.448 ± 0.030 - i0.026 ± 0.024) GeV 
is also found. In the bottom sector, one broad and one narrow poles are found at about 
(5.536±0.029-z0.117±0.043) GeV and (5.842±0.022-z0.018±0.010) GeV, respectively The 
coupling constants of the generated states to the relevant coupled channels arc calculated. 
They are consistent with the results in the pole structure analysis. In the (0, 1) channel, 
the width of the isospin violated decays D*i^{2317) L>+7r° and B;^(5729) 5°7r° are 
calculated. They are about 8.69 and 1.54 keV, respectively. Finally in the (|, 0) channel, 
the decay ratio ri/(ri -|- r2) for the higher narrow state is also estimated. 



Acknowledgments 

We are very grateful to M.F.M. Lutz, H.-Y. Cheng, J. A. Oiler and Y.-L. Shcn for valuable 
discussions. This work is partially supported by the NSFC grant Nos. 90103020, 10475089, 
10435080, 10447130, CAS Knowledge Innovation Key-Project grant No. KJCX2SWN02 and 



13 



Key Knowledge Innovation Project of IHEP, CAS (U529). 



[1] 
[2] 
[3] 
[4] 
[5 



[6 

[7] 
[8] 
[9 

[lo: 
[11 

[12 
[13 
[14 
[15 
[16 



BABAR Collaboration, B. Aubert et al., Phys. Rev. Lett. 90 (2003) 242001. 
CLEO collaboration, D. Besson et al., Phys. Rev. D 68 (2003) 032002. 
Belle Collaboration, P. Krokovny et al., Phys. Rev. Lett. 91 (2003) 262002. 
Behe Collaboration, Y. Mikami et al., Phys. Rev. Lett. 92 (2004) 012002. 



E. Vaandering, hep-ex/0406044|, in Proceedings for the XXXIXth Rencontres de Moriond 
(QCD and High Energy Hadronic Interactions), March 28th to April 4th, 2004, La Thuile, 
Italy 

BABAR Collaboration, B. Aubert, et al., |hep-e x70408067 , at the 32nd International Confer- 
ence on High-Energy Physics (ICHEP 04), Aug. 2004, Beijing, China. 
S. Godfrey, Phys. Lett. B 568 (2003) 254. 

W.A. Bardeen, E.J. Eichten, C.T. Hih, Phys. Rev. D 68 (2003) 054024. 

P.Colangelo,F. De Fazio, Phys. Lett. B 570 (2003) 180. 

Y.-B. Dai, C.-S. Huang, C. Liu, S.-L. Zhu, Phys. Rev. D 68 (2003) 114011. 

Fayyazuddin, Riazuddin. Phys. Rev. D 69 (2004) 114008. 

M.A. Nowak, M. Rho, I. Zahed, Acta Phys. Polon. B 35 (2004) 2377. 

T. Lee, I.W. Lee, D.P. Min, B.-Y. Park, hep-ph/0412"210| 

P.Colangelo,F. De Fazio, A. Ozpineci, Phys. Rev. D 72 (2005) 074004. 



X.-H. Guo, H.-W. Ke, X.-Q. Li, X. Liu, S.-M. Zhao, lhep-ph/0510146 



T. Matsuki, T. Morii, K. Sudoh, hep-ph/0510269, at the Xlth International Conference on 
Hadron Spectroscopy (HADRON05), August 2005, Rio de Janeiro, Brazil. 
[17] J. Lu, W.-Z. Deng, X.-L. Chen, S.-L. Zhu, hep-ph/0602167| 
[18] H.-Y. Cheng, W.-S. Hou, Phys. Lett. B 566 (2003) 193. 
[19] Y.-Q. Chen, X.-Q. Li, Phys. Rev. Lett. 93 (2004) 232001. 
[20] H. Kim, Y. Oh, Phys. Rev. D 72 (2005) 074012. 

[21] M. Nielsen, R.D. Matheus, F.S. Navarra, M.E. Bracco, A. Lozea, |hep-ph/0509 13H at the 
Workshop on Light-Cone QCD and Nonperturbative Hadron Physics 2005 (LC2005), July 
2005, Cairns, Australia. 

[22] K. Terasaki, |hep-ph/0512285| at the Workshop on Resonances in QCD, July 2005, Trento, 



14 



Italy. 

[23] Z.-G. Wang, S.-L. Wan, |hep^ph/0602080l 

[24] A. P. Szczepaniak, Phys. Lett. B 567 (2003) 23. 

[25] T. Barnes, F.E. Close, H.J. Lipkin, Phys. Rev. D 68 (2003) 054006. 

[26] E.E. Kolomeitsev, M.F.M. Lutz, Phys. Lett. B 582 (2004) 39. 

[27] J. Hofmann, M.F.M. Lutz, Nucl. Phys. A 733 (2004) 142. 

[28] Y.-J. Zhang, H.-C. Chiang, P.-N. Shen, B.-S. Zou, |hep^h/060427T| 

[29] E. van Beveren, G. Rupp, Phys. Rev. Lett. 91 (2003) 012003. 



[30] Z.-G. Wang, S.-L. Wan, hep-ph/0603007 



[31] T.E. Browder, S. Pakvasa, A.A. Petrov, Phys. Lett. B 578 (2004) 365. 
[32] J. Vijande, F. Fernandez, A. Valcarce, Phys. Rev. D 73 (2006) 034002. 
[33] P. Bicudo, |hep-ph/0512'04T| 



[34] W. Wei, P.-Z. Huang, S.-L. Zhu, Phys. Rev. D 73 (2006) 034004. 
[35] M. Nielsen, Phys. Lett. B 634 (2006) 35. 



[36] X. Liu, Y.-M. Yu, S.-M. Zhao, X.-Q. Li, ihep-ph/0601017 

[37] Behe Collaboration, K. Abe et al., Phys. Rev. D 69 (2004) 112002. 

[38] FOCUS Collaboration, J.M. Link, Phys. Lett. B 586 (2004) 11. 

[39] M.E. Bracco, A. Lozea, R.D. Matheus, F.S. Navarra, M. Nielsen,Phys. Lett. B 624 (2005) 217. 
[40] J.A. Oiler, E. Oset, Nucl. Phys. A 620 (1997) 438; (Erratum) ibid. 652 (1999) 407; 
[41] J.A. Oiler, E. Oset, J.R. Pelaez, Phys. Rev. Lett. 80 (1998) 3452; Phys. Rev. D 59 (1999) 
074001. 

[42] J.A. Oiler, E. Oset, Phys. Rev. D 60 (1999) 074023. 
[43] N. Kaiser, Eur. Phys. J. A 3 (1998) 307. 
[44] J. Nieves, E. Ruiz Arriola, Phys. Lett. B 455 (1999) 30. 
[45] V.E. Markushin, Eur. Phys. J. A 8 (2000) 389. 

[46] F.-K. Guo, R.-G. Ping, P.-N. Shen, H.-C. Chiang, B.-S. Zou, accepted for publication in Nucl. 
Phys. A Ihep-ph/OSOgPSO] . 



[47] J.A. Oiler, Phys. Lett. B 477 (2000) 187. 

[48] J.A. Oiler, Phys. Lett. B 500 (2001) 263. 

[49] J.A. Oiler, Nucl. Phys. A 725 (2003) 85. 

[50] L. Roca, E. Oset, J. Singh, Phys. Rev. D 72 (2005) 014002. 



15 



[51] G. Burdman, J.F. Donoghue, Phys. Lett. B 280 (1992) 287. 

[52] M.B. Wise, Phys. Rev. D 45 (1992) 2188. 

[53] T.-M. Yan, H.-Y. Cheng, C.-Y. Cheung, G.-L. Lin, Y.C. Lin, H.-L. Yu, Phys. Rev. D 46 (1992) 
1148. 

[54] J. Hofman, M.F.M. Lutz, Bucl. Phys. A 763 (2005) 90. 

[55] A.M. Badalyan, L.P. Kok, M.I. Pohkarpov, Yu.A. Simonov, Phys. Rept. 82 (1982) 31. 

[56] J.A. Oiler, Phys. Rev. D 71 (2005) 054030. 

[57] P. Cho, W.B. Wise, Phys. Rev. D 49 (1994) 6228. 

[58] R. Dashen, Phys. Rev. 183 (1969) 1245. 



